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Abstract 

We solve the one dimensional Feshbach-Villars equation for spin-1/2 particle subjected 
to a scalar smooth potential. The eight component wave function is given in terms of the 
hypergeometric functions and via a limiting procedure, the wave functions of the step 
potential are deduced. These wave functions are used to test the validity of the boundary 
conditions deduced from the Feshbach-Villars transformation. The creation of pairs is 
predicted from the boundary condition of the charge density. 
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The eight component relativistic wave equation for spin-1/2 particle, called the Feshbach- 
Villars equation (FV~), has been constructed and used to solve physical problems [H El 11] • 
The hydrogen atom is the first problem solved by Robson and Staudte [U [2]; they found the 
same bound-state energy as the Dirac equation but the wave functions are different. Recently, 
Robson and Sutanto have solved the Compton scattering problem and found that the cross 
section is given, like in Dirac theory, by the Klein-Nishina formula [5] . The same authors have 
also calculated the transition probabilities for the Balmer and Lyman a-lines of hydrogenic 
atoms and compared them to the Dirac and Schrodinger results [B]. 

Following the idea of a previous paper [7J, we study the one dimensional FV\ equation 
for a particle subjected to a step potential V(x). As we know the above problem is trivial 
in quantum mechanics but in relativistic quantum mechanics it emerges the famous Klein 
paradox. In addition to the known literature in this framework [5], we limit ourselves in this 
paper to illustrate the problem more clearly and justify the phenomenon of pair creation, and 
we left the problem of interpretation to the specialist researchers. The boundary conditions 
for the eight component wave function for the case of the step potentiel are unknown. Then, 
in order to bypass this problem, we take V(x) as a smooth potential [Sl HHl [II]- The analytic 
solution of the FV-1/2 equation with the smooth potential is given. We deduce, via a limiting 
procedure, the wave functions of the step potential and derive the transmission and reflection 
coefficients. The comparison with the Dirac coefficients constitutes the test to our calculations. 
The appropriate boundary conditions for the step potential are extracted from the Feshbach- 
Villars transformation. The boundary condition for the charge density is also evaluated. The 
validity of these boundary conditions is tested using the wave functions of the step potential. 
At the end, we discuss the boundary condition of the charge density and compare its predictions 
to those obtained from the transmission and reflection coefficients. 

The Dirac equation can be written in a second-order form as 

[( 7 ^) 2 + m 2 ]$ = 0, (1) 

where 7^ are the Dirac matrices, is the minimally coupled derivative and í>(^ 1; 9? 2? V3? ^4) 
is the four component wave function. This equation can be also written in the Klein-Gordon 
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form as follows 



m 2 )l 4 + V^]<3> 



0. 



(2) 



where a^ u = §[7^,7"] and = d^Ay — d^A^. The last term in Eq. (T5]) represents the spin 
interaction with the external electromagnetic field F^. By analogy with the equation of spin-0 
particle, the second order equation ([2]) is called the Klein-Gordon equation for spin-1/2 particle 
(KG\) [U E] . For spin-0 particles, this equation is reduced to a Klein-Gordon type equation. 

In order to linearize Eq. (j2J) to a first order equation in time, the Feshbach-Villars lineariza- 
tion procedure is used to transform the four component wave function $ to an eight component 
wave function ip [TJ [2]. The eight component wave function ip satisfies a Schrodinger type 
equation [2] 



with 



H 



(T 3 + IT 2 ) 



1 

2m 



d 



-D 2 l 4 + -a^ 



+ m(r 3 <g> 1 4 ) + eA l 8 , 



(3) 



(4) 



where r 2 , r 3 are the Pauli matrices, (g> is the Kroneker(direct) product, D = d + ieA is tridi- 
mensional minimally coupling and (A, A ) is the electromagnetic potential. 

The Hamiltonian H is pseudo-Hermitian H = r^H^ T4 and the inner product is 



where r 4 = r 3 <g> 7 . This inner product is indefinite in sign, it takes positive or negative values. 
Then, the (FV-1/2) wave function space is not a Hilbert space and its dimension is twice that 
of the Dirac solution space [2]. 

In the Weyl representation of the gamma matrices, Eq. ([3]) separates into two four- 
component equations with the Hamiltonians 



Ht = (r 3 + it 2 ) 



H„ 



(r 3 + ít 2 ) 



1 

2m 

1 
2m 



D 2 1 2 + iea.(E + iB) + m(r 3 ® 1 2 ) + eA l 4 , 



-D 2 1 2 - iea.(E - iB) 



m(Tz ® I2) + eA l 4 , 



(5) 



(6) 



where E, B are the intensities of the electromagnetic field and (t(t\, t 2 , r 3 ) are the Pauli matrices. 
The Hamiltonians H^ and H v satisfy also a Schrodinger type equation 



d 



(7) 



d 

where ip^ and ip v are four-component wave functions deflned by their components as 



(8) 



Ipf = (^1,^2,^3,^4 



1p V = (^5,^6,^7,^8 



and the eight component wave function is ip = (ip^,ip v ) • The Hamiltonians H% and H v , the 
wave functions ip^ and ip v transform to each other under spatial inversion, respectively [U [2] . 
In the Weyl representation of the gamma matrices, the density p is defined [2] as follows 

p = ipip where ip = V> f r 5 , r 5 = n <g> (r 3 (g) 1 2 ), 

where p and r 5 are given by Eqs. (62), (63) in Ref. [2]. 
We define the one dimensional current j by 
1 
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- dip dip 
^O— - -r^Oip 
ox ox 



—AipOip, 
m 



where O 



(t 3 + ít 2 ) 



u. 



The values of p and j are independent of the representation and they satisfy the continuity 
equation. p is interpreted as the charge density of the particle. The positive solution ip, 
the negative solution ip c , the charge density and the current are transformed by the charge 
conjugation as follows 



ip — > ip c = n ® 7°^ f , p 

In one dimension, Hamiltonians (JHI) , (J6 

V(x) are 



Pc = ~P 



and j 



Jc=3- 



Hr, 



(t 3 + ir 2 ) ® 



1 

2m 



d 2 



for spin-1/2 particle subjected to a scalar potential 



(r 3 + «r 2 ) <g> — 
2m 



dx 2 
dx'- 



zen 



-I2 + ien 



dx 

dV(x) 
dx 



+ m(T 3 ®l 2 ) +eV(x)U, (9) 



+ m(r 3 ®l 2 )+eV(:r)l4. (10) 



The terms ±zen 
potential. 



dx 



represent the interaction of the spin with the derivative of the scalar 



The stationary solution has the form ip(x,t) = e ip(x) for each component of the wave 



function and Eq. 
ential equations 
d 2 



of the Hamiltonian H^ is equivalent to the following four coupled differ- 



dx 2 



+ 2m 2 + 2meV(x) - 2mE 



ipi- 



d 2 ip 3 : dV(x) 



dx 2 



— le- 



dx 



(lp2 + Ipi) = 0, 



'ir 



d 2 



dx 2 

d^Pi 
dx 2 

d 2 ý 2 
dx 2 



+ 2m 2 + 2meV(x) - 2mE 
d 2 



"p2 - 



dx 



+ 



+ 



dx 2 
d 2 



dx 2 



+ 2m 2 - 2meV(x) + 2mE 
+ 2m 2 - 2meV(x) + 2mE 



4 . dV(x) , , , , 

dy(x) 



^3 - ie 
■04 — ie 



dx 



(02 + Ipi) 
(01 + V> 3 ) 



0. 



The difference of Eqs. (ÍTTT) - (ÍTÏÏI) and (ÍT2I) - fíTil) give, respectively. 

E - eV(x 



if>i ~ ^3 



rn 



E - eV(x) 



rn 



(4>i + fo), 
(^2 + ^4)- 



^2 - ^4 

Using these equations, the sum of Eqs. (fTT |) -( ÍT3j) and (fT2l)-(ÍT4l) give 
d 2 



dx 2 



d 



+ [E- eV(x)] z - m 1 



j-l+lE- eV(x)} 2 - m 2 



(V>i + ^3) +ie^^(ý 2 + ^) 

(ifo + + ze ^2r^i + ^3) 



0. 



The sum and the difference of the two last equations (ÍT71) . ( ÍTSl) give 



2 dV(x) 
m + ie- 



m 



le- 



dx 

d V(x) 
dx 



r e (x) = 0, 



VÁx) = 0, 



where 



(12) 
(13) 
(14) 

(15) 
(16) 

(17) 



(19) 
(20) 

(21) 



^l(x) =^+^2+^ + ^4, ip«(x) = + Ïp3 ~ ~ fa. 

We note that the differential equation (ÍTÏÏj) is the same as the Dirac equation gives (Eq. (207.7) 
in Ref. [T3j ) . If we take V(x) as a step potential V(x) = V)9(x), we have in Eqs. (ÍT9l) .(j20l) the 
potential V(x) and its derivative: V' (x) = dV ^ = V)5(x). In the configuration space, the delta 
Dirac potential has particular treatment of its boundary conditions [TD, [15] . In this paper, 
our goal is to find the wave functions without using boundary conditions. At the end, we use 
the obtained wave functions to test the validity of the boundary conditions deduced from the 
Feshbach-Villars transformation. For this, we use the scalar smooth potential V(x) defined as 



n*) = y(l+tanh|l) 



(22) 



where V and r are positive constants. In the limiting case r — > 0, V(x) — > V 9(x). It increases 
from the value V = for x = — oo to the value V = Vq for x = +00, the main rise occurring 
in the interval -2r < x < +2r : V(-2r) = 0.1192Vb, V (2r) = 0.8807V . 

In order to find the solution of the differential equations ÍÍT9|) , ff20l . we make the change of 
variable 

y = I(l-tanh|;), (23) 
which maps the interval x G ]— 00, +00 [ to y 6 ]0, 1[ . The new form of Eqs. (ÍTÏÏj) . 



arc 



1 



y 2 (i - y) 



, d 2 ^l(y) 

dy 2 



+ 



+ 



-y(l-y)(l-2y) 



dy 



.eV n 



(-S + eV y - eV ) 2 - m 2 + i— y(l - y) 

r 



r 6 (y) = 0, (24) 



■ 2n ,2 ^(1/) 1 n V1 9 ^Vf(y) 



dy 



+ 



eV 



(E + e y y - eV ) 2 - m 1 - i—y(l - y) 



^l(y) = 0. (25) 



The singularities of these differential equations are y = 0, 1, 00. Let us introduce the change 
V'l(y) = 2/"(l ~ yY f(y)i an d ^(y) = y"(l — yY9(y)i the last equations are reduced to the 
hypergeometric equation form 



1/(1 - 1/)^¥ + [(2^ + 1) - v&v + ^ + 2)] /y ' 



2/(1-2/) 



dtyv) 



[(2i/ + 1) - 3/(2i/ + 2/i + 2)] 



^(2/) 



1 1\2 U 

(/^ + ^+2) 



f///- ' ' Ji dy 

where i/ 2 = r 2 [m 2 — (E — eV ) 2 ] , /i 2 = r 2 (m 2 — £" 2 ), 
uf = 1 - 4r 2 e 2 V 2 - 4zreV . 
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(/Í + Z/ + -) 2 -- 2 - 
VP 2 ; 4 

„2„2t/2 



/(2/) = 0. (26) 



g(y) = 0, (27) 



1 - 4 r V V + AireVo and 



In comparison with the spin-0 case, the imaginary part in the expressions of v\ and v\ 
represent the effect of the spin. If we remove them, Eqs. (1261) . (1271) are reduced to the same 
equation as for the spin-0 case [7J. 

The general solutions of Eqs. (|26p . (|27|) are given in terms of the hypergeometric function 

+C 12 y~ v (\ - yY 2 F x (n - v + \ + n - v + - - 1 - 2v, y), 
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iPf(y) = D iy "(l-y)» 2 F 1 (// + i/+l + ^,/x + i/+Í-^,l + 2ï/,y) 

+L> 12 y-"(l - y)" 2 F!( M - i/ + - - -j, n - v + - + -j, 1 - 2u, y ), 

where C\, C 12 , í+ and _Di 2 are constants. We note that these solutions can be obtained directly 
from Eqs. ( 1241) . ( 125]) using a symbolic software [16j. We have chosen the parameters of the 
hypergeometric function in order to have an analogy with the spin-0 solution |7J . An equivalent 
solution with others parameters is given by Eq. (207.15) in Ref. [13] . 
We choose the regular solutions at the origin y = 

^l(y) = C, y v (l - yY 2 Fi(/i + v + 1 - |, /i + v + 1 + |, 1 + 2!/, y), (28) 

e d (y) = A y"(l - 2 F!(/x + */ + i + /i + + i - |, 1 + 2z/, y). (29) 

Then, the expressions of the components of the wave function tp^ can be deduced as follows: 
from the definition fl2T|) of ip^(y) and ipf(y), we have 

+ ^Áv) = 2(Vi + ^s), V>I(y) - ^ d (z/) = 2(^ 2 + 4 ), 



and using the Eqs. (fTBT) , (ÍTBl) we have 



^1,3(2/) 
^2,4(2/) 



X ± E-eV(y) 
m 

E-eV(y) 
m 



1 ± 



bPI(y)+^l(y)}, 
m(y)-^(y)], 



(30) 



where ip^(y) and 0|(y) are given by Eqs. (128]) , (1291) . the sign (+) corresponds to first index and 
the sign (— ) to the second index. 

The wave function ip v is calculated by the same method as ip^. We note that ip v can be 
also deduced from ip^ under spatial inversion. It satisfies Eq. ([8]) where the hamiltonian H v 
is given by Eq. (ÍTUl) in the one dimensional case. Its components satisfy the four coupled 
differential equations ífTT]) - (fT4"]) where the term (~3' e ~^ ) is replaced by the term (+^e^p) 
and the components (ipi,ip2, 03, ipi) by (05, ipQ, ipr, ip%, ), respectively. We give here the essential 
results: 

'E-eV(x) 



05 ~ ifa 
^6 - 08 




(31) 
(32) 



— + [£ - eV(a;)] 2 - m 2 - ie— — 

UiJU H>Jb 

d 2 r „ T ,, N1 2 2 aVíx) 

- + [E-e^)] 2 - m 2 + ,e^U 



^(x)=0, 
<(z) = 0, 



(33) 
(34) 



where 



ýfa) = ^5+^6 + ^7 + ^8, ^(^=^5 + ^7-^6-^8- (35) 

We note that ^ and i/j^ satisfy the same differential equations as ipf and respectively. Then, 
the solutions of the differential equations ( 1331) and ( 134"1) are similar to the solutions of Eqs. ( 1201 ) 
and (ÍTÏÏj) . respectively 

^(y) = C 2 ^(l-y)^ 2 F 1 ( /U + z/ + i + |, /U + z/ + i-|,l + 2z/, 2/ ), (36) 

<(</) = ^2 y"(l - yY 2*1(0 + v + ~ - y , // + z/ + ~ + y , 1 + 2u, y), (37) 

where C 2 and .D 2 are constants. 

Then, the expressions of the components of ip v can be deduced as follows: from the defini- 
tion (1351) of ipí(y) and ipt(y), we have 

%(v) + <(y) = 2(^5 + ^7), ^(y) - <(y) = 2(^6 + V> 8 ), 



and using Eqs. (lïïTj) . (1ÏÏ21) we have 

Aáv) 



x ± E-eV(y) 
m 

l ± E-eV(y) 



[r v (y) + Vr,(v)l 
my)-tá(v)l 



(3* 



where ^(y) and V'n(y) are gi ven by Eqs. (1ÏÏB1) . (13T|) . 

We study now the asymptotic behavior of the wave function when x — > ±oo. First, when 
x — ► —00 or y — ► 1, we have (1 — y) ~ exp(a;/r); we use the property of the hypergeometric 
function which links the y and (1 — y) argument, 

2^1(0,6,0,3/) = A 2^(^,6,0 + 6-0+1, l-í/) + J B(l-í/) c - a - b 2*1(0-0,0-6,0-0-6+!,!-;?/), 



with 



.4 



r(c)r(c-a-6) 
r(c-a)r(c-6) 



r(c)r(a + 6-c) 
TRr(6l 
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The corresponding constants of the waves functions é^ and éf are : 



A, 



r(2u+i)r(-2fi) 



T(2u+l)r(-2fi) 

r(^- M +i-f )r{u-»+í+% )' 



5, 



r(2i/+i)r(2/x) 
r( M +v+i-^)r( M +v+i+^) ' 



(39) 



T(2u+1)T(2ij,) 

r( M +^+i+f )r( M+i ,+i-f )' 



lim^^l, lim(l - yY = e^ x/r , lim(l - y)'» = e"^ /r , and 3^1(0, 6, c, 0) = 1. 

2/->l y->l j/-»l 

Thus, when 2 — ► —00 or y — > 1, the waves functions and ^| have the following behavior: 

éi(x) — ► A s e» x/r + B s e~^ /r , é?(x) — ► A d e^ x/r + B d e-" x/r . (40) 

; X — *— 00 s x^— 00 

Setting /i = —irki, with k\ = E 2 — m 2 , where fci is real positive 

éi(x) — ► A s e~ iklX + B s e iklX , éf(x) — ► A d e~ iklX + B d é klX . (41) 



x— ►— 00 



x— >— 00 



For the limit when x +00, or y — ► 0, 2 Fi(a, c, 0) = 1, \imy u = e ra/r , and lim(l —yY = 1- 

y^O y^O 

Then, the waves functions have the following behavior: 



-ux/r 



x— >+oc 



s x— >+oo 



-ux/r 



(42) 



Setting v = —irk 2 , with k\= (E — eVo) 2 — m 2 while k 2 is real for E < eVo—m or E > eVo+m 
and k 2 is imaginary for eVo — m < E < eVo + m. Then, the waves functions are 



r e (x) 



x— >+oo 



^ x , éï(x) • i 



ik^x 



x— >+oo 



(43) 



By the same method the asymptotic behavior of the waves functions é s and é^ are 

é s Jx) — > A d e- Ífeia: + 5 d e ifcl *, <(x) — > A s e~ lklX + B s e iklX , (44) 

' X— > — 00 ' X— >— 00 



V>"(a;) — ► e ífc2X , ^(z) — ► e ik2X . (45) 

We note here that for the plane waves and the other wave function which are not square- 
integrable the renrmaliztion condition takes the form [2] (ipE, Í>e') = iS(E — E'), while for the 
square-integrable wave function it takes the form (é,é) = ±1. 

The reflection and transmission coefficients can be calculated from the current density in 
the one dimensional case 

1 



J 



2im 



Tn dé dé 
ox ox 



A = 0, 



(46) 
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where O = 1 2 ® (73 + ir 2 ) ® 12- Using the last definition of the current and the incident wave, 
we find that the incident current is 

ki 

jinc = — [B*B d + B s B* d }. 
m 

The reflected current is evaluated using the reflected wave 

Jref = -— [A* s A d + A s A* d ] . 
TXX 



Then, the reflection coefficient is 



\jref \ _ \A* s A d + A s A d \ 



__ \jrej 1 _ ____ 1 ____ 

\jinc\ \B* s B d + B s B* d \ 
The transmission coefficient is evaluated in terms of the transmitted wave, 

jtr = — (h + k_) expi(k 2 - k\)x. 

If k 2 is real, 



. _ 2k 2 

Jtr — j 

m 

and the transmission coefficient T is 

rr, \jtr\ 2h 1 



(4Í 



\ Jmc \ h \B*B d + B s B* d \ 
If k 2 is imaginary, 

j tr = 0, T = and R=l, (49) 

in this case we have a total reflection. 

We consider now the limiting case when the smooth potential tends to step potential, i.e., 
when the parameter r tends to 0. 

For the first region x < 0, the limits of the coefficients (|39|) when r — > + are 

hm A s = kl ~ h - eV \ limB s =^M, (50) 
r^o+ 3 2k_ r^o+ 2k x y 1 

A h-k^ + eVo kx + k 2 -eV_ 
hm A d = , hm B d = , 51 

r^o+ d 2ki ' r^o+ d 2k_ y 1 



10 



and the waves functions and ipf(x) are 

ki — kz — eV 



= Ci0(-ar) ■ 
= D\9(-x) 



2k\ 

k\ - k 2 + eV 



exp(—ik\x) + 



2ki 



k\ + k 2 + eV 
2k\ 

k\ + k 2 - eV 



2k\ 



exp(+ik\x) | , (52) 
exp(+ik\x) > . (53) 



exp(— ik\x) + 

The terms +eV in the above relations are the contribution of the spin and if we remove them, 
Eqs. ( I52"l) .(l53"l) are reduced to the same equation as for the spin-0 case \7\. 
For the second region x > 0, the waves functions are similar to Eqs. (l4"3l) 



^(x) = C\9(x)e tk2X , i>l(x) = D\9(x)e 



D ik2X 



Then, the waves functions can be written in compact form for the two regions 



(54) 



^|(x) = C\9(-x) ■ 
ipf(x) = D\9(-x) 



k\ - k 2 eV 
2k\ ~2~k\ 



k\ - k 2 eVo 
2k\ 2k\ 



exp(—ik\x) + 
exp(—ik\x) + 



k\ + k 2 eV 
2k\ + 2~k[_ 

k\ + k 2 eV 
2k\ 2k\ 



exp(+ik\x)j+C\9(x)e ik2X , 

(55) 

exp(+ik\x)^+D\9(x)e ik2X . 

(56) 



Using the same method, the waves functions ip^ and ip^ in the two regions are 



= C 2 9(-x) ■ 
tf(x) = D 2 9(-x) 



k\-k 2 eV 



2k\ 2k\_ 

k\ - k 2 eV 
2k\ 2k~\ 



exp(—ik\x) + 
exp(—ik\x) + 



k\ + k 2 eV 



2k\ 2k\ 

k\ + k 2 eVo 
2k\ 2k\ 



exp(+ik\x)j+C 2 9(x)e lk2X , 

(57) 

exp(+ik\x)X + D 2 9(x)e ik2X . 

(58) 

In the case of the step potential, we note that the presence of the delta Dirac in the differential 
equations fTT9j) . (f20|) . (T33Ï) . (T34l ) implies that the functions (J55l) - (l58l) are continuous at x = and 
their derivatives discontinuous. Then, the final expression of the wave function ip(x) of the 
step potential is deduced using the expressions of the eight components (1301) . (1381) and the last 
relations (155D-(155D: 



ip(x) 



ip 2A (x) 

Í>5j(x) 
\ A,s( x ) ) 



1 ± 



E - eV 9(x) 



rn 



1 i>ï(y) + i>ï(y) ^ 

ip^(y) - ipf(y) 
r v (y) + ip d v (y) 
{ ^(y)-^(y) ) 



(59) 



n 



where the sign (+) corresponds to the first index and the sign (— ) to the second index and 

h - k 2 



6{-x) 



(Ci±D0 



2k x 



+ 



r v {y)±ip d v {y) =0{-x) 



(m _l n \ kl ~ , eV 0(ry n 



+ 



(n _l n + e ^0^ n \ 



exp(— ifcix) 
exp(+iA;ia;)| + 0(rr)(C a ± D x )é k2X ', 

exp(-ikix) 
exp(+ik lX ) \ + 6(x)(C 2 ± D 2 )e ik2X . 



We note also that the wave function of spin-0 particle [3, [T7] can be deduced from the wave 
function (159|) of spin-1/2 particle if we take only the two components (^i,^) an d remove the 
effect of the spin-1/2 term (eVó). This analogy with spin-0 particle does not exist for the Dirac 
wave functions. 

At the end, from Eqs. (]4Tl). (H8l) and (Í501) . (Í5TT) we deduce the reflection coefficient R and 
the transmission coefficient T for the step potential 
• For k 2 real positive (k 2 > 0) and E > eV + m, we have 



R 



(h - k 2 f - (eV o y 
(fci + k 2 f 



ihk 2 



(eVnÝ' 



(h + k 2 ? - (eV ) 



2 : 



and R + T = 1. 



(60) 



For k 2 real negative (k 2 < 0) and m < E < eVo — m, we have 
(h + hf - (eV ) 2 „ Ahk 2 



R 



T 



and R — T 



(61) 



(h-hf- (eV f ' (h-hf- (eV f 

which is the well known Klein Paradox. 

The coefhcients T and R given in the above relations coincide exactly with the Dirac 
ones ^ïïl [TB| [T5] . If we remove the effect of the spin-1/2 term (eVo) from the above relations 
( I60p . (161j) we find the reflection and the transmission coefncients of the spin-0 case [7]. 

Like in the spin-0 case [7], in the following we are going to look for boundary conditions 
using the Feshbach-Villars transformation. In the case of the step potential, the presence of 
the delta Dirac in Eq. (j2J) implies that the KG-1/2 wave function $ is continuous at x = and 
its derivative discontinuous. 
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The Feshbach-Villars transformation is defined for the eight components CcLSG cLS follows: 



V»i,(ar,í) 



^i(a;,í 
ip 2 (x,t 
ip 3 (x,t 
\ ýi(x,t 

( ip 5 (x,t 
ipe(x,t 
ip 7 (x,t 
\ ip 8 (x,t 



1 

V2 



(pi + 

<P% + 
<Pl ~ 



Í(| + ^)^i 



(62) 



1 

V2 



<Pa + 

P>3 ~ 
V ^4 - 



m^dt 



ÍeV)(f4 

, H + ieV)ip 4 J 



" 9 + f e y)^ 



(63) 



which can be written in the abridged two-component form 



( <p, ï 




V V>i+2 / 





+ á(á + ieVr )^ 
^-á(l+ íeF )^ 



(64) 



where the couple (j, k) takes the following values: (j, k) = {(1, 1), (2, 2), (5, 3), (6, 4)}, respec- 
tively, (ipj,ipj + 2) T are the eight components of the wave function ip with j = 1,2,5,6 and 
<Pk = <Pk(x, t) are the four components of the wave function $ with k = 1, 2, 3, 4. 
From the last definition (|64|) it follows that 

1 



Pk 



V2 



(iPj + ý 



J+2), 







^dt ~ eV ^ k 



m 
V2 



(ipj - ip j+2 ) . 



(65) 



(66) 



The stationary KG-1/2 wave function $ has the form $(x, t) = e lEt $(x) for each component 
and the last equation is written as 



m 



(67) 



(E - eV)p k = -j= (ipj - ip j+2 ) . 

The continuity at x = of the component p k of the KG-1/2 wave function defined in Eq. 
gives 

^•(0+) + iP J+2 (0 + ) = iP 3 (0-) + ip á+2 {0-). (68) 

The continuity at x = of the component p> k of the KG-1/2 wave function defined in Eq. (jS 
gives 



^■(0 + )-^- +2 (0 + ) 



E-eV 
E 
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ip á (o~)-ip á+2 (o-) 



(69) 



From Eqs. (j680 and (IBTIl) we can write the boundary conditions in the matrix form 



/ 



/ 



__ 

2E 



2E 



/ 



. __ 

\ 



2_B 



V ^j +2 (°~ 



(70) 



where j = 1, 2, 5, 6. 

For the following particular cases the boundary conditions are more simple: 



For E = ^: 



^■(o + )=^ +2 (o- 



^■ +2 (o + ) = ^(o-). 



This case will be discussed at the end. 
For E = eV : 
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^•(o + ) = ^ +2 (o + ) = -[Vi(o-) +^ +2 (o-)]. 

The boundary conditions ( 170|) are between the components of the wave function ip. It can 
be written for all the eight components as follows: let us introduce the following notations 

= = + + ips + if>a, 

3=1,2,5,6 

fal = ^'+2 = ^3 + ^4 + ?p7 + ^8, 

j=l,2,5,6 

from relation (1681) . the sum of the eight components of the wave function satisfies 



(71) 



^/(o + ) + ^//(o + ) = ^j(o-) + V//(o-). 

From relation ( 1691) . the difference between -0/ and ipu satisfies 

E-eV 



M® + ) ~ M(0 n 



E 



<M0") - <Mo- 



(72) 



From relations íjTTl) and (1721) we can also write the above boundary conditions in the matrix 
form 

<mo- 



i _ £_. sYsl 

2E 2E 



eVp i _ __ 

2E 2E 



V/(0 + ) 
^//(0 + ) 

Using the boundary conditions (1701 . we find that the charge density p is discontinuous 

E — eV 



(73) 



P(0 H 



E 



-p(0- 



(74) 
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it gives the charge sign of the transmitted particle as a function of the energy E, the potential 
Vo, the charge e and the charge density of the incident particle. We note that the same relation 
( I74p can be also obtained for spin-0 particle. In the case of the Dirac equation, the density is 
continuous at x — 0. Then, the Dirac density is a probability density and cannot be a charge 
density . 

At the end, we verify that the wave functions ( l59l) of the step potential satisfy the boundary 
conditions (1701) . (1731) . (1741) and the current(H6l) is continuous at x = 0. 

The boundary conditions (17T1) . (1721) can be also interpreted by analogy with electromagnetic 
waves (when they traverse two different regions) as follows: the sum of the components of the 
wave function ip s = ipi + ipn is continuous like the tangential component of the electric field but 
the difference of the two components ipd = ipi — ipn is discontinuous like the normal component 
of the magnetic field. 

Now we discuss the boundary condition ( 1741) of the charge density and compare its predic- 
tions to those obtained from the transmission and reflection coefficients for the three following 
cases(we assume for all cases [IH] that eV > 2m and E > m i.e. k\ > 0): 

1. E>eV + m 

is real and < E ~^ Vo < 1 ; from Eq. (1741) the charges densities p(0 + ) and p(0~) have the 
same sign, i.e. the charge of the transmitted particle has the same sign as the charge of the 
incident one. For k^ > 0, the transmission and the reflection coefficients are given by relation 
( 160|) . For E » eV + m: p(0 + ) ~ p(0~), T ~ 1 and R ~ 0. The two methods give the same 
results. 

2. eV - m < E < eV + m 

k<i is imaginary and from Eq. (174|) we consider two cases: 

• For eV < E < eV + m: we have E ~^ Vo > 0; the charges densities p(0 + ) and p(0~) have the 
same sign. 

• For eVo — m < E < eV : we have E ~^ v ° < 0; the charge densities p(0 + ) and p(0~) have 
opposite sign, this means that we have creation of particle-antiparticle pairs near the step 
barrier if the potential is strong enough. 
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On the other hand, from Eq. (1491) we have a total reflection R = 1 and T = and the 
wave function (loT?]) is decreasing (evanescent wave) in the second region P3H [20] . The only case 
for which we have a total reflection from Eq. (1741) is for E = ëV§. p(0 + ) = 0. 

3. m < E < eV — m 

k 2 is real and E ~^ Vo < 0, from Eq. (Í74"I) the charges densities p(0 + ) and p(0~) have opposite 
sign. This means that we have also creation of pairs near the step barrier. For k 2 < 0, the 
transmission and the reflection coeflicients are given by Eqs. ( 161 p . Then, the creation of pairs 
in the Klein Paradox [T5] is proved from the boundary condition of the charge density (1741 . We 
note that Guang-jiong Ni et al [2ÏÏ1 EI] have discussed the Klein paradox, for the spin-0 case, 
using the current and the charge density in the two regions but they haven't used the boundary 
condition of the charge density. 

Let us now study an interesting particular case for the spin-0 particle[7J: 
• For the particular value E = from the interval m < E < cVq — m (we have assumed 
that eVo > 2m): we have k 2 = ±&i, T = 1 and R = i.e. the incident particle is transmit- 
ted to the second region (called the resonance transmission) but from relation (I74p we have 
p(0 + ) = — p(0 - ): this means that we have creation of pairs near the step barrier and the trans- 
mitted particle is the antiparticle of the incident one. The result of the case k 2 = — ki need an 
interpretaion. 

In summary, in order to find the wave functions of the step potential without the use 
of boundary conditions, we introduce the smooth potential as an intermediate stage. Then, 
we solve the one dimensional Feshbach-Villars equation for spin-1/2 particle subjected to the 
smooth potential. The eight-component wave function is given in terms of the hypergeometric 
functions. In the limiting case r — > 0, the wave functions of the step potential are deduced 
in each region. The transmission and reflection coefficients are identical to the Dirac ones. We 
have also an analogy between the wave functions and the transmission and reflection coefficients 
of the spin-1/2 and the spin-0 particles. Boundary conditions relative to the step potential are 
extracted using the Feshbach-Villars transformation and the continuity of the KG-1/2 wave 
function at x = 0. The main result is that boundary conditions for the step potential are: 
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• the sum of the eight components ij) s — ipi + ipn is continuous: 

^(o + ) = ^(o-), 

• the difference of the two components ipd 

Mo + ) - 

• the charge density p is discontinuous: 

P(0 + ) = 

and for E ~^ v<) < we have creation of particle-antiparticle pairs if the potential is strong 
enough. Then, the number of particles becomes variable and this implies that we must use 
quantum field theory [Hl [13 122] - 

At the end, we note that we have omitted the singular solution in order to make comparison 
with the Dirac results. We propose to study the contributions of this solution separately to 
this paper. 
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